Spin polarized STM spectra of Dirac Fermions on the surface of a topological insulator 
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We provide a theory for the tunnehng conductance G{V) of Dirac Fermions on the surface of 
a topological insulator as measured by a spin-polarized scanning tunneling microscope tip for low 
bias voltages V . We show that GiV) exhibits an unconventional dependence on the direction 
of magnetization of the tip and can be used to measure the magnitude of the local out-of-plane 
spin orientation of the Dirac Fermions on the surface. We also demonstrate that if the in-plane 
rotational symmetry on the surface of the topological insulator is broken by an external field, then 
G{V) acquires a dependence on the azimuthal angle of the magnetization of the tip. We explain 
the role of the Dirac Fermions in this unconventional behavior and suggest experiments to test our 
theory. 

PACS numbers: 71.10.Pm, 73.20.-r, 72.25.-b 



I. INTRODUCTION 

Topological insulators in both two and three dimen- 
sions (2D and 3D) have attracted a lot theoretical and 
experimental attention in recent years^"'^,. It has been 
shown in Refs. 4- 6 that such 3D insulators can be com- 
pletely characterized by four integers vq and J^i,2,3- The 
former specifies the class of topological insulators to be 
strong (i^o = 1) or weak (y^ — 0) while the latter in- 
tegers characterize the time-reversal invariant momenta 
of the system given by M = i^2^2, i^3^3)/2, where 

^1,2,3 are the reciprocal lattice vectors. The latter inte- 
gers 1^1,2,3 are basis dependent while the former, vq, is 
a basis-independent invariant for a given material. The 
weak topological insulators (WTI) are adiabatically con- 
nected to Anderson insulators whereas the strong topo- 
logical insulators (STI) are not; consequently the topo- 
logical features of STI are robust against the presence 
of time-reversal invariant perturbations such as disor- 
der or lattice imperfections. STI exhibit a host of novel 
phenomenon such as induction of magnetic monopoles 
in the presence of an electric charge near its surface^, 
presence of topologically protected Fermion modes inside 
dislocations^, and the possibility of realization of mag- 
netic control over electrical conduction in its junctioni2. 
It has been theoretically predicted^ and later experi- 
mentally verified^ that the surface of a STI has an odd 
number of Dirac cones whose positions are determined 
by projection of M to the surface Brillouin zone. The 
position and number of these cones depend on both the 
nature of the surface concerned and the integers J^i.2,3- 
For several compounds specific surfaces can be found 
which will host a single Dirac cone near the F point of 
the 2D Brillouin zone^iii. Such a Dirac cone at the sur- 
face of a topological insulator is described by the Dirac 
Hamiltoniani^ 
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where (?(/) denote the Pauli (identity) matrices in spin 
space, ip{k) = {ip^{k),ipi{k))'^ is the electron annihilation 
operator, k — (kx,ky) is the 2D wave vector, z denote 
the unit vector normal to the topological insulator sur- 
face, and fi is the chemical potential. The properties of 
these surface Dirac electrons has been studied in detail 
in recent years. In particular, these Fermions are ex- 
pected to exhibit spin-momentum locking which predicts 
their spin to be along the surface. This phenomenon has 
been extensively studied by recent spin-resolved ARPES 
studiesiiii^. However, in many cases, such studies lead 
to a finite, albeit small, probability of the measured spin 
polarization to point along the z direction^^. It is un- 
clear at the moment whether this contribution comes 
due to possible additional higher order terms in the Dirac 
Hamiltonian of the surface electrons which can spoil the 
spin-momentum locking property, or due to the presence 
of additional scattering potentials induced by point-like 
or step-like defects. Such spin- momentum locking for the 
surface electrons has also been studied indirectly through 
scanning tunneling microscope (STM) studies near step 
edgesiS. However, the surface of these insulators has not 
been studied using spin-polarized STMi^. 

In this work we develop a theory for the tunneling cur- 
rent of the Dirac Fermions on the surface of a topological 
insulator as measured by a spin-polarized STM tip. We 
derive explicit expressions for the tunneling current I{V) 
and the tunnehng conductance G{V) — dI{V)/dV and 
show that these quantities exhibit an unconventional de- 
pendence on the direction of the tip magnetization. Such 
an unconventional nature of the tunneling current and 
conductance originates from the fact that Fermions on 
the surfaces of topological insulators, in contrast to con- 
ventional electrons obeying the Schrodinger equation, ex- 
hibit spin-momentum locking which does not allow a free 
choice of their spin quantization axis (which is fixed along 
the z direction as evident from Eq. ([1])). We demonstrate 
that the analysis of G{V) provides a direct method of 
measurement of the local out-of-plane spin component 
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FIG. 1: (Color online) Schematic representation of the pro- 
posed experimental setup. The polar (a) and the azimuthal 
(/3) angles of the tip magnetization Mo are shown. The ap- 
plied electric and magnetic field are along z and x respectively. 



of the Dirac Fermions on the surface. Further, we show 
that if the in-plane rotational symmetry on the surface 
of the topological insulator is broken by an external field, 
then G{V) acquires a dependence on the azimuthal an- 
gle of the tip magnetization. We substantiate our theory 
and demonstrate these unconventional features by a com- 
putation of G{V) for Dirac Fermions in the presence of 
a crossed electric (in-plane along x) and magnetic field 
(out of plane along z) as shown in Fig. We suggest 
experiments which can test our theory. 

The organization of the rest of the paper is as follows. 
In Sec. ini we obtain the general expressions for I{V) 
and G{V) as measured by a spin-polarized STM tip at 
low bias voltages V. This is followed by Sec. IIIII where 
we provide an explicit calculation of the tunneling con- 
ductance G{V) for Dirac Fermions on the surface of a 
topological insulator in the presence of a crossed in-plane 
electric and out-of-plane magnetic field. We discuss pos- 
sible experiments and conclude in Sec. IIVI 



II. TUNNELING CURRENT DUE TO A 
SPIN-POLARIZED STM 



The experimental situation that we intend to describe 
is schematically shown in Fig. ([Ij. The STM tip is placed 
atop the surface of the topological insulator and has a 
magnetization Mq as shown in Fig. ([T|). To understand 
the unconventional behavior of the tunneling current in 
such a system, let us briefly review the current obtained 
by such a spin-polarized tip when placed atop the surface 
of a conventional magnetic material with a magnetization 
mo(r). In this case, the STM current is given by the well 
known expressior 
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where / is the Fermi distribution function, E*^ and E^ de- 
note the energy levels of the tip and the surface states re- 
spectively, and M^i, denotes the matrix element for over- 
lap between the tip and the STM wave functions which 
is determined by the Bardeen tunneling formulai^S 

2m Jg ^ 

-C(r,2RV;,*(r,^)], (3) 

where r denotes only the (x, y) coordinates. Here V'Jj and 
11)1 denote the wave functions of the electrons in the tip 
and on the sample surface respectively, z is taken to be 
the direction perpendicular to the surface, and the inte- 
gral is to be performed over a plane which lies between 
the STM tip and the sample surface and is parallel to the 
sample surface. 

For a conventional magnetic sample, one can always 
choose the local spin quantization axis to be along the 
direction of magnetization of the tip. Consequently the 
wave function of the tip electrons can be written as 



^^(r;z) = (l,0)x^(r;z), 



(4) 



where the precise form of X/i (r) is difficult to obtain the- 
oretically and it is usually modeled to have the same or- 
bital symmetry of the electron for the atom at the STM 
tip. For the rest of this work, we shall assume an s-wave 
symmetric wave function for the tip electrons; all the 
results obtained can be extended to other orbital sym- 
metries in a straightforward manner using the methods 
of Ref. 19. Note that such a choice of the spin quanti- 
zation axis of the tip is possible since both the surface 
and the tip electrons obey Schrodinger equations which, 
in contrast to the Dirac equation, do not specify a choice 
of the spin quantization axis. For the surface electrons, 
the wave function can be written as 

rAv,z) = (cos(^,(r)/2),sin(0,(r)/2)e'*''W)e.(r;z) 

(5) 

where the angles 9^(y,z) and 4'i,{r,z) are polar and az- 
imuthal angles with respect to the tip magnetization, and 
(r) denotes the spatial part of the surface electron wave 
function whose specific form depends on the systemi^i^i. 
Using Eqs. (O, dH) and ([5]), a straightforward analysis 
following Ref. |2l| yields 



Co cos(6l^(rt)/2)^^(rt). 



(6) 



where rt = (rt,Z() denotes the position of the tip, and 
the constant cq depends on the details of the tip wave 
function whose expression for a spherically symmetric s- 
wave tip is given fa yiSi^OiZi 



Co 



(7) 
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where A'' is the normalization constant, n is the decay 
length of the tip wavefunction, and m is the mass of the 
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tip electrons. Substituting Eqs. ([6]) in Eq. ([2]) one obtains 
a more familiar expression for the current 

I{V) = Io\co\^pt{ps{eV;rt) + ms{eV,rt)), (8) 

where V is the applied voltage, Iq = 2ne/h, pt is the 
density of states for the tip electron which is assumed 
to be constant within the range of experimental applied 
voltages, and the integrated local density of states ps = 
Ps{eV; rt) and magnetization nis = ms{eV; rt) are given 
by 



J doj [ficj) 
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-^)\Urt)\' 




-fiuJ + eV)] 




-t^)|C.(rt)pcos(0,(rt)). (9) 



Note that if we restrict ourselves to low bias voltages and 
low temperatures for which the contribution to comes 
from a narrow range of v for which the relative angle of 
magnetization between the sample and the STM tip is 
approximately a constant, we obtain 

liV) = Io\cofptPs{eV; n) [1 + cos(0(rt))] (10) 

which is the well-known expression for the current for a 
magnetic material as measured by a fully polarized STM 
tip. 

The central result of our work, which we now derive, 
is that Eq. ([T0| needs to be significantly modified for 
the calculation of I{V) for a topological insulator. For a 
topological insulator, the electrons on the surface are de- 
scribed by the Dirac Hamiltonian given by Eq. ([T|) with 
some possible additional warping terms, impurity poten- 
tials, and/or external electric and magnetic fields. In all 
such cases, the spin quantization axis of the electrons is 
fixed to be along z. In this situation, for the system of 
polarized STM tip near the surface, once we choose to 
describe the Dirac Fermions by Eq. ([T]) (with possible 
additional terms), it is no longer possible to fix the spin 
quantization axis of the system along the magnetization 
of the STM tip. Thus the tip wave function has to be 
represented as 

i>lir;z) = (cos(a^(r,z)/2),sin(a^(r,z)/2)e^'^-('--^)) 

^xA^;z), (11) 

where Q;^(r, z) and /3^(r, z) denotes the polar (as mea- 
sured with respect to z) and azimuthal angles of the 
magnetization (as measured with respect to an in-plane 
rotational symmetry breaking field along x which will be 
specified later) as shown in Fig. ([Ij. Note that in the ab- 
sence of an in-plane symmetry breaking field, as shown in 
Fig. (2)), the azimuthal angle can be chosen to have any 



arbitrary value. The wave function of the Dirac electrons 
on the surface of a topological insulator in real space can 
be written as 

£\r;z) = (^Ti(r;z),^Ti(r;z)) 

= {l,K,{r,z)e-''^'^^'^^^^)^Jl{r;z) (12) 

where ip'^l{r;z) and ^pjl{r;z) represent the spin- up and 
spin-down component of the Dirac electrons and we have 
used V'Ji(r; z)/'0^J(r; z) = K^(r, z) exp[-z77i,(r, z)]. We 
will provide explicit expressions for these wave functions 
for a specific case in the next section. 

Substituting Eqs. ^ and ^ in Eq. ([3]), one then 
obtain the matrix element Af^^ entering the expression 
of the current I{V) in Eq. @ to be 



M^i, = co-0t^(r;z) cos(a^(r, z)/2) + ^^.(r, z) 
X sin(a,,(r, 2)/2)e'['3M(r,^)-').(r^^)] 



(13) 



which yields, assuming constant density of states for the 
tip and energy independence of the tip magnetization 
angles a and /?, 

I{V) = /olcolVt PdieV;rt) + Pz{eV;rt)cos{a) 



PdieV;rt) 



Pz{eV;rt) 



Pm(eV",/3;rt) 



+/9,„(ey, /3;rt)sin(Q;) 
^ fcL; [fiu) - fiuo + eV)] 



dLO [f{uj) - f{LO + eV)] 

X^<5(£;.-..)|V'^T_I(^^)|2[i_^2(^^)] 

j doj [ficj) -f{oj + eV)] Y S{E, - ^) 
x|Vt.T'^-(rt)cos[/3-77,(rt)] (14) 



Using Eq. [141 the tunneling conductance Giy") at T = 
can be obtained from Eq. (fT4| as 



G = Gq\cq\'^ pt\p'd{(^V ■,rt) + p'^{eV;Yt)cofi{a) 
+pL(e^>/3;rt)sin(a) 

V 

X cos[,9 - 77,.(rt)], (15) 
where Go = ^-Kf? jfi. 



Pd(eT^;rO 
P™(eV',/3;rt) 
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Eqs. (|T4)) and ([15]) represent the central result of this 
work. First, we note that the tunneling conductance dis- 
plays an unconventional dependence on both the polar 
and azimuthal angles of the tip magnetization. In par- 
ticular, G{V) may exhibit a dependence on the azimuthal 
angle /3 of the magnetization of the tip which never oc- 
curs for STM spectra of conventional magnetic samples. 
Second, we note that we do not expect this unconven- 
tional dependence to show up when the tip is placed on 
top of a pristine topological insulator surface without any 
defects and/or in the absence of external fields. In this 
case, since the spins of the Dirac electrons point along 
the plane due to spin-momentum locking, one expects 
Ki,{rt) = 1 which leads to pz — 0. Further, since the 
sample is completely isotropic, the sum over all the sur- 
face states with a given energy E^, in Eq. ((T?| will make 
the current vanish. This can be seen from the fact that 
such a sum must be independent of the orientation of 
the X and y axis on the surface of the TI; therefore shift- 
ing rjiy rj^ + c should not change any of the results. 
Hence, the quantity pm in Eq. must be indepen- 
dent of the choice of (3 and must therefore vanish. This 
leads to I{V) = Io\co\'^ptPd{eV;rt) which is independent 
of both a and /3. Finally, we would like to point out 
that the presence of a non-zero component of Sz due to 
warpingi^ or other effects such as presence of defects^^ 
renders K,y(rt) ^ 1 and hence leads to a cos(q!) depen- 
dence of I{V) and G{V). Thus measuring G{V) with a 
spin-polarized STM provides a direct way of measuring 
the real-space out-of-plane spin polarization of the Dirac 
electrons on the surface of a topological insulator. We 
shall discuss this point in more details in Sec. IIVI 

In the next section, we compute G{V) for the Dirac 
electrons in the presence of a crossed electric (along x) 
and magnetic field (along z) and show that G{V) displays 
all the unconventional features discussed above. 



III. CROSSED ELECTRIC AND MAGNETIC 
FIELD 

In this section, we compute G{V) for the geometry 
shown in Fig. ([T}, i.e., in the presence of a constant 
magnetic field B along z and an electric field E along 
X. Throughout this section, we shall assume that the ap- 
plied magnetic field is weak enough so as not to orient the 
magnetization of the tip along B. For a stronger mag- 
netic field, one needs to apply the magnetic field along the 
direction of the magnetization of the tip. Note that for 
the topological insulator, this just amounts to changing 
B ^ Bz = B cos(a) since the in-plane component of the 
field appears to be a gauge shift for the Dirac electrons 
and does not have any infiuence on their properties-'^^. So 
all of the analysis of the present section will hold in this 
case, with B — )• Bcos{a). 

We first look for the solution of the wave function 
ipj^^^{r,z) — V'J-^_^(r) exp(— K^z), where Kz is the decay 
length of the wave function in real space which depends 



on the work function of the topological insulator surface. 
For estimating G{V), we need the value of the wave func- 
tion at the position of the tip, i.e., at z — zt. Thus it 
is possible to absorb exp(— K^zt) in the definition of Go 
and one does not need to know the precise value of Kz 
to estimate G. The solution for ?/'J-f (r) and the cor- 
responding energies E^ for the present configuration is 
well-knowrilfli2.2. 'We will assume that the electric field 
is smaller than the magnetic field [E < vpB); then one 
can use a Lorentz transformation to go to a frame where 
the electric field is zero, solve the Landau level problem in 
that frame, and then Lorentz transform back to the origi- 
nal frame. The energy eigenstates are labeled by the Lan- 
dau level index n and the transverse momenta ky and can 
be expressed in terms of A = arccosh[(l — i?^/t>|,i3^)~^/^] 
for 71 ^ by 

<i,t« = [cosh(A/2)x„-i(e)+isinh(A/2)x„(e)]e''^''^ 
= [cosh(A/2)x„(^)-*sinh(A/2)x„-i(0]e"^''^ 
E^^k, = Sgn(n)SoCosh(A)"3/2 

X Vl^l + CzB cosh(A) - hkyE/B, (16) 

where Eq = HvfIb', Xn = exp(-eV2)g„ (g)/ ^2 1"! |n|!, 
Hn denotes Hermite polynomials, Ib — \/ hc/eB is the 
magnetic length, Cz — 9^ P% I {^v'j^e.) , g is the gyromag- 
netic ratio, /is is the Bohr magneton, and ^ is given by 

^^{x-kyll)/^+^\E/{vFB). (17) 
From Eq. (|16p , one obtains the parameters k and rj to be 

_ ^cosh^(A/2)x^(g) + sinh^(A/2)x?,I^ 

y'cosh2(A/2)xLi(0 + sinh2(A/2)x2 (0 ' 
rjn.kyix) = arctan[tanh(A/2)x„(C)/x«-i(C)] 

-|-arctan[tanh(A/2)x„-i(e)/Xn(0] ■ (18) 

The corresponding expressions for the n = Landau level 
are given by 

En^ky = -gusB/ cosh{X) - hkyE/B, 

<fc„t(r) - e*'=='^*sinh(A/2)xo(0, 

= e^'=='^cosh(A/2)xo(e), (19) 

which leads to kq — coth(A/2) and 770 — — 7r/2 indepen- 
dent of ky . We note here that for typical topological insu- 
lators gpBB/Eo ^ 10""* for _B ~ IT and so the Zeeman 
term can be safely neglected in all numerical estimates 
of G. 

To obtain G, we substitute Eqs. ([T5|) and (|18p in 
Eq. dTSl) . In what follows, we choose the STM tip to 
be directly on top of the origin of coordinates so that 
rt — (0,0). We note that corresponding to any given ap- 
plied voltage V, the solution of En,k = eV constitutes, 
in principle, an infinite number of pairs (n,ky). How- 
ever, the Gaussian factor exp(— ^2/2) in the expression 
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FIG. 2: (Color online) Plot of G/Go as a function of the ap- 
plied voltage eV/Eo for a = (top left panel) and a — tt (top 
right panel) for E/vpB = 0.1 and |coppt/So = 1. The bot- 
tom panels show the variation of the wave functions |Xi-i (?) 1^ 
(bottom left) and |xn(OI^ (bottom right) corresponding to 
the solution {n,ky) of eV = Eniky) which contributes maxi- 
mally to G. See text for details. 




E/(v^B) E/(v^B) 

FIG. 3: (Color online) Plot of G/Go as a function of the 
applied electric field E/vpB for several applied voltages 
eV/Eo — 1.4 (red solid line), 1.8 (green dotted line), and 
2.4 (blue dashed line) for which the maximally contributing 
{n,ky) pair corresponds to n = 1, 2 and 3. The left (right) 
panel corresponds to a = (tt). We have taken \co\'^pt/Eo — 1 
in all the plots. See text for details. 



of ^(r) and the discrete values of n ensures that out 
of all these pairs, the one with the minimal value of \ky\ 
provides the most significant contribution to G. This is 
particularly valid for small E <^ vpB which is also the 
regime we are interested in. In this regime, the tunneling 
conductance G{V) as a function of V , as computed using 
Eqs. ([H]), dUl) and ([15]), is shown in the top left (right) 
panel of Fig. ^ for E/vpB = 0.1 and a — (tt). From 
Eq. (fT5)) . we find that for both a = and a = tt, pj„ = 0, 
so that one has 

G{V;a = 0) = Go\co\^ PtY.S{E, - eV^^lm^ , 

G(y-a = i^) = Gn\c^\^ptY,5{E,-eV)\rPllm^. 

(20) 

The nature of G(T^) can now be understood from Eqs. 
(dni and (Uni). From these equations, we find that as V 
is changed so that it deviates from the energy of a Lan- 
dau level, the lowest value of ky = ky which is a solution 
to eV — En,ky starts to increase. This leads, through 
a change of £,{k^)^ to a non-monotonic change in the 
wave functions 1-0^3(0)^ (for a = 0) and |i/'JJ(0)p (for 
a — -k) since these wave functions involve a product of the 
Gaussian factor exp(— ^^/2) and the Hermite polynomi- 
als 7?„_i(0 and Hn\0- Note that for E/vpE = 0.1, the 
predominant contribution to |'0^J(O)P[IV'J!(O)I^] comes 
from the coefficients of the cosh(A/2) terms and thus in- 
volve Hn-i{£,) [Hn{S)] for a = [tt]. In fact, the plots 
of Xn-i[^{ky)] {Xn[^{ky)]) as a function of V in the lower 
left(right) panel of Fig. ^ (for a = (tt)) show that the 
non-monotonic behavior of G{V) can be understood in 
terms of the variation of these wave functions with the 



change of ky as a function of V . We note that the spec- 
tral features in GiV) are expected to change drastically 
with the direction of the tip magnetization as is evident 
from comparing the left and right panels of Fig. ^ . 

Next, we fix the applied voltage eV/Ea and study the 
behavior of G'(y) as a function of the electric field E/vpB 
for several values of eV/Eo for a = (left panel) and tt 
(right panel) as shown in Fig. ([3]). The nature of these 
plots can again be understood from the variation of the 
wave functions tpJl{Q) and ^'^^(O). For example, let us 
consider the plot corresponding to cV/Eq — 1A (red solid 
lines in Fig. ([3])) for which the maximally contributing 
{n,ky) pair corresponds to n = 1. As one increases 
E/vpB, we find from Eq. PT)) that the value of ^ de- 
creases since the decrease in ky overcompensates for the 
increase in E leading to a net decrease in ^ in the regime 
where E/vpB ^ 1. Consequently, for a — 0, where 
G depends on \^Ji,{0)\^ - lxo(OP - exp(-e2), we find 
an increase in G{V). For a — n, where G depends on 

\^JM? - 1X1(01' - \exp{-e/2)H^{0?, a similar in- 
crease in f leads to a non-monotonic behavior of G. The 
plots of G for other values of eV/Eo shown in Fig. ([3]) 
can be understood in a similar manner. 

Finally, we consider the variation of the tunneling con- 
ductance as a function of the polar (a) and azimuthal (/3) 
angles corresponding to the direction of the tip magne- 
tization for an applied voltage cV/Eq = 1.4 and electric 
field E/vpB = 0.1 as shown in Fig. (g]). We note that 
G{V) shows maximal variation with /? for a = 7r/2 as 
expected from the expression of pj„ in Eq. (ITS|) . Such 
a variation is highlighted in Fig. ([5]) where the plot of 
G'/Go (where G' = dG/dfi) is shown as a function of /3 
for a = 7r/2. We note that since G'(/3 = 0,a = 7r/2) ~ 
sin(?7(rt; T^)), the intercept of this plot is a measure of 
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FIG. 4: (Color online) Plot of G/Go as a function of the angles 
a and /3 for E/vfB = 0.1, \co\^pt/Eo = 1, and eV/Eo = 1.4. 



file of the surface Dirac electrons near surface defects 
such as point impurities and step edges. The other ex- 
periment we suggest involves measurement of the STM 
spectra with a tip whose magnetization points along the 
surface of the topological insulator (a = 7r/2) in the pres- 
ence of a crossed electric and magnetic field. Our central 

G 
Co 




the local relative phase between the spin-up and spin- 
down components of the surface Dirac Fermions. Such 
a variation of the tunneling conductance with the az- 
imuthal angle of the tip magnetization is a novel feature 
of the topological insulators which originates from fixing 
the spin quantization axis of the surface Dirac electrons 
along z, and it has no analog in conventional magnets. 
We note that the measurement of such a variation in an 
experiment does not require a change in the direction 
of the tip magnetization. It can simply be achieved by 
changing the direction of the electric field E which breaks 
rotational symmetry since /? is always defined with re- 
spect to the direction of E. 



IV. EXPERIMENTS 

The experimental test of our work involves measure- 
ment of the STM spectra of a topological insulator using 
a magnetic tip. We note that such measurements can 
be used to map the spatial out-of-plane spin-orientation 
profile of the surface Dirac electrons. In particular, we 
suggest measurement of the spectra with the STM tip 
magnetization pointing along z {a — Q) and — z [a = tt). 
The local z component of the spin of the Dirac electron 
can then be obtained from 



FIG. 5: (Color online) Plot of G'(V)/Go (where G' = dG/dp) 



G{V-a^Q)-G{V]a^TT) 
G{V\a = Q) + G{V;a^TT) 



(21) 



as a function of /3 for E/vfB • 
and eV/Eo = 1.4. 



0.1, \corpt/Eo = 1, a = 7r/2. 



We note that such a measurement would provide direct 
information about the real-space out-of-plane spin pro- 



prediction is that the tunneling conductance will display 
an oscillatory behavior with the change of the azimuthal 
angle /3 (which can be varied by changing the direction 
of the applied electric field for a fixed tip magnetization) 
as can be seen from plot of dG/d/3 as a function of /3 in 
Fig. (O. We also note that the value of G"(/3 = 0) pro- 
vides direct information about the local relative phase of 
between the spin-up and spin-down components of the 
electron wavefunction on the topological insulator sur- 
face. 

In conclusion, we have presented a general theory for 
the STM spectra of Dirac electrons on the surface of a 
topological insulator as measured by a magnetic STM tip 
and have shown that such a spectrum has unconventional 
features not commonly seen in its counterpart in conven- 
tional magnets. We have identified the reason for such 
an unconventional spectrum to be the fixation of the spin 
quantization axis of the surface Dirac electrons along the 
z direction. We have also explicitly computed the STM 
spectra for the surface Dirac electrons in the presence of 
a crossed electric and magnetic field which demonstrates 
this unconventional behavior. 
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